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We investigate the phenomenology of top squarks at the Large Hadron Collider (LHC) in a 
supersymmetric model where lepton number is identified with an approximate U(1)r symmetry 
in such a way that one of the left chiral sneutrinos can acquire a large vacuum expectation value 
(vev) and can play the role of the down-type Higgs. This A-symmetry allows a subset of trilinear 
A-parity violating interactions, which determine the collider phenomenology of this model in a 
significant way. The gauginos are Dirac particles and gluinos are relatively heavy in this class of 
models. The model contains a right handed neutrino superfield, which gives a tree level mass to one 
of the active neutrinos. An order one neutrino Yukawa coupling also helps enhance the Higgs boson 
mass at the tree level and results in a very light bino-like neutralino 0^) with mass around a few 
hundred MeV, which is a carrier of missing (transverse) energy (JSr )• The model can accommodate 
two rather light top squarks, compatible with the observed mass of the Higgs boson. The lighter 
top squark (ti) can decay into tx®, and thus the signal would be similar to the signal of top quark 
pair production at the LHC. In addition, fully visible decays such as t 2 —> be + can give rise to 
interesting final states. Such signals at the LHC combined with other features like a heavy gluino 
could provide a strong evidence for this kind of a model. Our analysis shows that 575 (750) 

GeV and mj 2 ^ 1.2 (1.4) TeV can be probed with 5<r statistical significance at the 13 TeV LHC 
with 300 (3000) fb _1 of integrated luminosity. Finally, we observe that in the presence of super-light 
carriers of JUt, the so-called ‘stealth’ top squark scenario may naturally appear in our model. 

PACS numbers: 14.80.Ly, 12.60.Jv, 14.65.Ha, 11.30.Pb 


I. INTRODUCTION 

The discovery of a Higgs boson at the Large Hadron Collider (LHC) with a mass around 125 GeV 0 , 0 | is of 
immense importance in high energy physics and, in particular, in the context of electroweak symmetry breaking. 
However, in spite of its enormous success over the years, the Standard Model (SM) of particle physics suffers from 
several drawbacks. From a theoretical perspective, the naturalness problem remains a troublesome issue in the 
framework of SM. Supersymmetry (SUSY) renders an elegant solution to this problem and has become the most 
popular choice for physics beyond the standard model (BSM) to date. Nevertheless, searches for superpartners by 
the LHC collaborations (ATLAS and CMS) in pp collisions at the centre-of-mass energies of = 7 TeV and 8 
TeV have shown no significant excess 03 over the SM background. This has put stringent lower limits on the 
superpartner masses in many different SUSY scenarios. Recent experimental analyses within the framework of a 
simplified phenomenological minimal supersymmetric standard model (pMSSM) have set a lower bound of 1.7 TeV 0] 
for comparable masses of the gluino and the first two generation squarks. 

On top of that, finding a Higgs boson with a mass ~ 125 GeV and the non-observation of any signals of physics 
beyond the SM have severely constrained many supersymmetric scenarios which are otherwise very well motivated. 
In view of this, models with Dirac gauginos and U(1)r symmetry have become popular as they can significantly 
lower the current exclusion bounds on the first and second generation squarks and at the same time can address the 
125 GeV Higgs boson even in the presence of lighter top squarks. Other virtues of such scenarios include significant 
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suppression of flavor- as well as CP-violating; effects. Assorted versions of models with Dirac gauginos and U(\)r 
symmetry can be found in the literature |6 h67| . 

On the other hand, experiments in the neutrino sector have firmly established the fact that neutrinos have tiny 
masses and non-trivial mixings. Non-vanishing neutrino masses and mixings [GSUtiI are very important indications 
of new physics. An interesting question to investigate is whether models with Dirac gauginos could also provide 
explanation for observed neutrino masses and mixings. As we shall describe later on, the introduction of a right- 
handed neutrino superfield with an appropriate R-charge and Yukawa coupling ‘/’ can give rise to a small neutrino mass 
of the right order at the tree level. At the same time, an order one ‘/’ generates an additional tree level contribution 
to the Higgs boson mass. Thus, physics in the Higgs sector and the physics in the neutrino sector become intimately 
related in this model. This gives an opportunity to look into the Higgs sector through the neutrino-window and 
vice-versa. The additional tree level contribution to the Higgs boson mass also opens up the possibility of having 
rather light top squarks in the spectrum. This can also ameliorate the situation with so called ‘naturalness’ that 
is somewhat compromised in popular SUSY frameworks like the MSSM. There, the top squark masses are pushed 
to higher values (~ O(TeV)) to ensure a Higgs boson as heavy as observed at the LHC. Another very interesting 
outcome of this scenario is to have a very light bino-like neutralino, also identified as the lightest supersymmetric 
particle (LSP), with a mass in the range of a few hundred MeVs. The scenario violates R-parity. Thus, such a light 
neutralino LSP could decay to SM fermions. However, the decay length turns out to be much larger than the collider 
dimension [72]]. Hence it would essentially contribute to missing transverse energy (MET; 0t)- 



The R-charges are identified with the lepton numbers in such a way that the left chiral sneutrino vevs can be large, 
and are not constrained by the Majorana mass of the neutrinos. Thus, the sneutrino can play the role of a down 
type Higgs field. There also exists a subset of R-parity violating operators, mixings between the neutrinos and the 
neutralinos, as well as, between the charged leptons and the charginos. Once U (l)fi symmetry is invoked, the gauginos 
cease to have Majorana masses. However, they can acquire Dirac masses which requires additional chiral superfields 
living in the adjoint representation of the SM gauge group. It is somewhat crucial in the context of the present work 
to note that R-symmetric models also prohibit the traditional trilinear scalar couplings (‘A' terms) and the Higgsino 
mass parameter (/r term). To generate a fi term, one needs to incorporate two more chiral superfields R u and Rd 
with appropriate R-charges. 

The main motivation for having a right handed neutrino superfield is to have a tree level neutrino mass [52J. 
However, such a simple extension has enormous implication for the lightest CP even Higgs boson mass [52], [63?] and 
for the dark matter sector [58 ]. 

The squarks in the present context carry a non-zero R-charge (R = 1) and hence a non-zero lepton number since 
lepton number is identified with the R-charges. Top squarks can naturally be light in this model and can have novel 
signatures at the LHC. In the present work, we take the obviously natural direction of connecting to the top squark 
sector which is very much in the focus of the current LHC programme and thus could be put to test in a straight¬ 
forward way. This work presents for the first time the collider implications of the very characteristic top squark sector 
of the scenario under consideration. Various possibilities in the decays of both top squarks are discussed in detail. A 
novel final state in the form of bbe + e~ is highlighted where the final state objects can, in principle, be reconstructed 
to the mass of the heavier top squark. The decay of the top squark to be + differs from typical R-parity violating 
MSSM decay modes, and is typical for this model. This also provides us with an interesting handle, using which the 
reach of t 2 can be enhanced significantly, so much so, that the enhanced rate could also lead to its discovery even 
before its_hghter cousin. 

As for 1 1 , pair production of top squarks and their subsequent decays might lead to signal similar to top quark pair 
production and provide important information on the model and, in particular, on the scenario with an order one 
neutrino Yukawa coupling ‘/’. A characteristic difference in the signal is in the form of a somewhat softer when 
compared to similar mj values in the MSSM. This is because of the presence of a MeV neutralino LSP, which is again 
a salient feature of our scenario. It is also demonstrated how various decay modes of t\ remain simultaneously open 
thus necessitating a thorough analysis of the experimental data from the 13 TeV run of the LHC. 

The paper is organised as follows. A brief description of the model is presented in section [TT] The neutral scalar 
sector of the model and its characteristic features are discussed in section mu Section m describes the electroweak 
gauginos in general. We start with a generic discussion of the neutralino sector and its role in generating the tree level 
neutrino mass both in U(1)r conserving case as well as taking into account mild violation of this R-symmetry. We 
discuss the possibility of mixings among the neutralinos and the neutrinos. Later on, we discuss the chargino sector 
of this model and the corresponding mixing between the charginos and the electron. The focus area of this work, that 
is to say, the top squark sector is described in section El Expressions for the decay rates in various relevant modes are 
presented. The latest bounds on the the masses of the top squarks as reported by the LHC collaborations are also 
discussed. The model is incorporated in SARAH (v4.4.1) [73W75I ]. In section ED we present a few benchmark points 
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that reflect the characteristic decay patterns of the two top squarks and are found to be instrumental in shaping the 
interesting signatures at the LHC. Section lAHl is devoted to the actual simulation study of the signals and the most 
relevant backgrounds using event generators. Estimations of the reaches in the masses of the the top squarks are 
also presented. In section IVIIII we briefly analyse the issue of the ‘stealth’ top squark which arises naturally in our 
scenario. We summarise with some concluding remarks in section HXl 


II. THE C/(l) fl -LEPTON NUMBER MODEL 

We minimally extend an Absymmetric model, first discussed in EM [46[ , by a single right handed neutrino superfield 
N c jhl]. Along with the chiral superfields of the minimal supersymmetric standard model (MSSM) superfields, H u , 
Hd, , Qi, Uf, Df, Li, Ef, the model contains two ‘inert’ doublet superfields R u and Rd with opposite hypercharges. 
To prohibit spontaneous J?-breaking and hence the emergence of Abaxions, the scalar components of R u and Rd are 
barred from receiving any nonzero vev. This is why R u and Rd are labeled as ‘inert’. Similarly, the scalar component 
of N c does not acquire any nonzero vev. The Dirac gaugino masses can be constructed with the introduction of chiral 
superfields, living in the adjoint representation of the SM gauge group. A singlet S is needed to form a Dirac mass for 
the U( 1) gaugino, a triplet T under SU(2)l is required to have a Dirac mass for the SU(2)l gauginos and similarly an 
octet O under SU(3)c must be there to generate the Dirac gluino mass. The U(1)r charges of the chiral superfields 
along with their SM gauge quantum numbers are shown in table [I] 


Superfields 

SU(3) C , SU(2) l , U(1 )y 

U(1)r 

Q 

(3, 2, |) 

1 

Ui 

(3, 1, -|) 

1 

m 

(3, 1, 1) 

1 

Li 

(1, 2, -1) 

0 

Ef 

(1, 1, 2) 

2 

Hu 

(1. 2, 1) 

0 

Hd 

(1, 2, -1) 

0 

Ru 

(1, 2, 1) 

2 

Rd 

(1, 2, -1) 

2 

S 

(1, 1, 0) 

0 

T 

(1, 3, 0) 

0 

6 

(8, 1, 0) 

0 

N c 

(1, 1, 0) 

2 


TABLE I: Chiral superfields with the SM gauge quantum numbers and U(1)r charge assignments. 


Note that the scalar components transform in the same manner as their respective chiral superfields whereas the 
fermions have Abcharge one less than that of the corresponding chiral superfields. Following reference [45[, we also 
identify the lepton numbers of the component fields to the negative of their Abcharges. Such an identification leaves the 
lepton number assignments of the SM fermions unchanged from the usual ones while the superpartners acquire non¬ 
standard lepton numbers. As mentioned in the introduction, it is quite conspicuous that the left chiral sneutrino vevs 
can be large since they do not become constrained by the lepton number violating Majorana neutrino masses |45j. As 
a, result, the sneutrino can play the role of a down type Higgs field. It is now possible to integrate out the superfields 
R u and Hd, which simplifies the superpotential and the scalar potential considerably. At this point all the three 
sneutrinos can acquire substantial vevs. However, without any loss of generality, one can always choose a basis in 
which only one of the sneutrinos get a non-zero vev, which we choose to be the electron type sneutrino 1 , whereas, the 
vevs of the other two sneutrino fields are zero. Thus, the electron sneutrino ( v a , a = 1(e)) plays the role of a down 
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It should be noted at this point that the formulation, though independent of this kind of a choice of a particular basis, could have crucial 
impacts on the actual observables at the experiments. We justify our choice later in this paragraph. 
















4 


type Higgs field. With this basis choice and the assumptions of [52| {63J, the superpotential takes the following form; 

W = + v u H u R d + fL a H u N c + \ s SH u R d + 2\ T H u TR d 

- M r N c S + W', (1) 

W = £ AYY + £ 23 kL'. 2 L' 3 E' k c 

6=2,3 fc=l,2,3 fc=l,2,3 

+ £ KjkKQ'A 0 ’ ( 2 ) 

j,fc=l,2,3;6=2,3 

where y^j is the up-type Yukawa coupling, y u is the Higgsino mass parameter consistent with the i?-charge assignments. 
‘/’ represents the neutrino Yukawa coupling, A s and A t are the trilinear couplings and finally Mr is the coefficient 
of the bilinear term N C S. 

Note that, for simplicity in this work we have neglected the terms nN c SS, r]N c , £ N c Tr{TT ) and ( N c Tr(00 ) from 
the superpotential. As long as y ~ Af| ugY , k, ~ 1 and vev-s of the scalar components of S and T are very small 
(~ 0(1CU 4 ) GeV, as we shall consider later), we do not expect any significant changes in the analysis and the results 
presented in this work. 

The prime in W' indicates the mass basis for the down type quarks and charged leptons. When the electron 
sneutrino gets a vev, the first two terms in W' give masses to the down type charged leptons and quarks. However, 
SU (2) invariance ensures that the LLE operator cannot generate a mass for the electron. The electron mass can be 
generated from higher dimensional operators suppressed by a mass scale as discussed in [45] . Such operators would, in 
principle, contribute to the masses of other charged leptons as well. However, these contributions will be subdominant 
compared to the contribution to their masses from the LLE operators. This makes it a natural choice for the electron 
sneutrino to acquire a non-zero vev. The other two terms in W' include all the trilinear R-parity violating but lepton 
number preserving terms in this model. 

In a realistic supersymmetric model, one needs to incorporate soft SUSY breaking terms such as the gaugino and 
the scalar mass terms. The Lagrangian consisting of the Dirac gaugino mass terms [45| can be written as 

Cugino = J d 2 e^[V 2/si W la S + 2V2k 2 tr (W 2a T) + 2sft k 3 tr(W 3Q 6)] + h.c., 

(3) 

where W' a = \ a + 6 a D' is a spurion superfield parametrising D-type SUSY breaking. Wi a ’s are the field strength 
tensors containing the gauginos of the MSSM vector superfields. The D-term vev generates Dirac gaugino masses 
which can be schematically written as Mf 1 = m <D A > , where sy's are the order one coefficients and A is the scale of 
SUSY mediation. 

Similarly, .R-conserving but soft SUSY breaking terms in the scalar sector can be generated from a spurion superfield 
X, where X = x + 9 2 Fx [45J. The non-zero vev of the R-term generates the scalar soft terms. In the rotated basis 
where only the electron type sneutrino acquires a vev, the soft SUSY breaking terms are given by 

Vsoft = m 2 Hu HtH u + ml d R] l R d + rnl a LlL a + Y m\L\L b + M 2 n N^N c 

6=2,3 

+ m\V R J, Ri + m^QjQi + m lY Ri URi + m 2 gd\idRi + m 2 s S^S 

+ 2m Y tr (T^T) + 2r?ZQtr(0^0) — ( ByLH u L a + h.c.) + (tsS + h.c.) 

+ i& s (S 2 + h.c.) + Mtr(TT)+h.c.) + R 0 (tr(00)+h.c.). (4) 

It is important to note that the scalar singlet tadpole term ( tsS ) is suppressed [42j in the scenarios with Dirac 
gaugino masses and that is what we will consider in the present context. With this short description of the theoretical 
framework we now proceed to describe the scalar and the fermionic sectors of the model in appropriate details. 


III. THE NEUTRAL SCALAR SECTOR AND THE STANDARD MODEL-LIKE HIGGS BOSON 

In this section we discuss the CP-even scalar sector, followed by a rather important discussion on the lightest 
CP-even mass eigenstate. The scalar potential receives contributions from the F-term, the D-term, the soft-SUSY 
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breaking terms and the dominant quartic terms generated at one loop and can be written down as 


V = Vf + Vb + Koft + bone—loop ■ 


( 5 ) 


From the scalar potential and the subsequent minimization equations, one can now write down the CP-even scalar mass 
matrix in the basis (h Rl vr, Sr, Tr), where the subscript R indicates the real parts of the corresponding superfields. 
Both Rd and N c carry //-charges of two units and hence gets decoupled from the CP-even scalar mass matrix. In the 
//-symmetric scenario the elements of CP-even 4x4 scalar mass matrix are given by j52j 


(Mf)u 

(M 2 S ) 12 

(Ms) 13 

(M 2 s)l4 

(Ms) 22 

(Ms) 23 

(Ms) 24 

( M s) 33 
(Af|)34 

(Ms)** 


(9 +g ) v 2 gin 2 p + (■ fM R v s ~ Bnl)(tan/3) 1 + 2SX u v 2 sin 2 /3, 

f 2 v 2 sin 2/3 + Bg a L - ^ 9 + 9 ^ sin 2/3 - fM R v s , 

2A| vs'V sin /3 + 2g u Xsv sin /3 + 2XsXtvvt sin /3 -f \[2g M 3 v sin /3 — fM R v cos /3, 
2X r vtv sin /3 + 2g u X rv sin /3 + 2XsXrVsv sin /3 — V2gM 2 v sin /3, 

( a 2 + n' 2 ) o 

-———-c 2 cos 2 /3 + (fM R vs — tan/3 + 2SX v v 2 cos 2 /3, 


-V2g'M°v cos /3 — fM R v sin /3, 
y/2gM 2 v cos /3, 

v 2 sin 2 /3 XsXtVtv 2 sin 2 /3 

l^u^S 

XsXtv 2 sin 2 /3, 

-li u X T — sin 2 /3 - XsX T vs — sin 2 /3 
Vt vt 


ts g'M^v 2 cos 2/3 fM R v 2 sin 2/3 
\/2ws 2ns 


gMf_v^_ 

\[2 vt 


cos 2/3. 


( 6 ) 


The S ’s appearing only in (Mj)n, (Mj) 12 and (Mj )22 quantify the dominant one-loop radiative corrections 2 to the 
quartic potential coming from the terms |5A u (|fI M | 2 ) 2 , i<5A 3 1| 2 1| 2 and ^5X V (|zy ,| 2 ) 2 where 


S X u 
SX 3 
SX U 


3 yf 

I6n 2 
5 A|, 
32t r 2 

M 

167T 2 


In 


In 


In 


fm.^m 

'*2 


5 A^ 

m( 

m^\ 

V TO 2 


) + 

167T 2 

v 2 ) 

(m 2 T \ 


X 2 

1 n 

to| 

\ + 

\ v 2 ) 


32t r 2 V 

u 2 




5A£ 

!n( 


V 2 

V TO 2 


) + 

167T 2 

. v 2 ) 


\4 
A S 

167T 2 



\4 

As In 


107T 2 


(S) 


(7) 


ms,rriT are the singlet and triplet soft masses while the singlet and the triplet vevs are denoted by vs and vt, 
respectively [22j]. g',g are the U(l)y and SU{2)l gauge coupling constants, M^,M 2 are the Dirac bino and wino 
masses, respectively, tan/3 = v u jv a , where v u is the vev of the up-type neutral Higgs field and v a represents the vev 
of the electron type sneutrino. The ellipsis at the end of each expression stands for missing subdominant terms. 

In this work, we study a simplified scenario in which the singlet and the triplet vevs are very small. These effectively 
decouple the corresponding scalar fields. Thus, the CP-even scalar mass-squared matrix turns out to be a 2 x 2 one 
and can be written down in a compact form as 

M 2 x = Mf sin 2 /3 + £ cot /?, 

M 2 2 = —£ + ^Mf(a — 1) sin 2/3 = M 21 , 

M 22 = £ tan /3 + Mf cos 2 /3, (8) 

where a = 2 ^ and £ = fM R vs — B^l- As long as M\ > Mf, where M\ = j s the CP-odd Higgs 

mass, we find that the tree level upper bound on the lightest CP-even Higgs boson mass is [ 52 } 

to); ^ [Mf cos 2 2/3 + f 2 v 2 sin 2 2/3] . (9) 


2 See reference [S3 , l45ll for a detailed discussion. 























6 


Clearly, this result is very interesting since for a large neutrino Yukawa coupling, / ~ 0(1), the Higgs boson mass 
receives a large tree level enhancement. This additional tree level contribution Tree = / 2 u 2 sin 2 2/3 grows at 

low tan/3 and becomes significant for order one neutrino Yukawa coupling [5211631]. The resulting enhancement 
could play a significant role in lifting the Higgs boson mass to 125 GeV. Furthermore, this additional contribution 
ameliorates the ‘naturalness’ (pertaining to the mass of the Higgs boson) issue in the MSSM. However, this tree level 
contribution gets diluted at large values of tan /3. There, the one loop quartic corrections EH can come into play 
and can substantially enhance the Higgs boson mass in the presence of order one couplings, A s and At, as shown 
in equation [7) Thus, even for larger values of tan/3, one can easily find a Higgs boson as heavy as observed at the 
LHC experiments, when the top squarks are relatively light. Figure Q] illustrates the region in the plane of and 



200 400 600 800 1000 1200 1400 1600 

m q(GeV) 


FIG. 1: Allowed region in (toj-- mj-) plane compatible with 124.7 GeV < rrih < 126.2 GeV Higgs mass after taking 
into consideration the full one loop corrections. The grey bands indicate the values of and my 2 ruled out by the 
LHC (see text for more details). The narrow vertical strip over 197 GeV < < 205 GeV refers to the ‘stealth’ top 

squark regime. 


mj 2 compatible with 124.7 GeV < mn < 126.2 GeV and various slices of A s over the ranges shown. This takes 
into account the one-loop corrections computed in the effective potential approach to the neutral scalar potential as 
implemented in SARAH (v4.4.1) [73l - f75| . For this scattered plot we use tan/3 = 23, M = M 2 D = 1.2 TeV, M 3 D = 
1.5 TeV, n u = 200 GeV, / = 1, v s = v T = 1CT 4 GeV, B(jl l = -(200 GeV) 2 , t s = (174 GeV ) 3 and vary A s in the 
range 0.8 < As < 1.3. We also vary the soft scalar masses (mg) 33 , (m 2 )33 in the range —8 x 10 6 (GeV ) 2 to 8 x 10 6 
(GeV ) 2 to vary the top squark mass. We observe that higher values of the superpotential coupling A 5 (hence larger 
At, as At = Astan$vv (®w being the usual weak mixing angle), chosen to fit the neutrino mass at the tree level) 
provide larger corrections to the Higgs boson mass at the loop level and hence the requirement of having multi-TeV 
top squarks can be avoided. The largest values considered for and mj 2 in this plot are close to 1660 GeV and 
1830 GeV, respectively. 

As can be seen from this figure, relatively light top squarks with mj — mj 2 may be generic to our scenario (though 
such a situation could attract more aggressive constraints from the LHC). Furthermore, both of them can have sub- 
TeV masses simultaneously and can still be consistent with the observed value of the SM-like Higgs mass if As > 1. In 
the absence of any appreciable chiral mixing (left-right ( L-R) mixing), top squarks as heavy as values considered here, 
could only raise the Higgs boson mass up to around 114 GeV in the MSSM. Hence the model under consideration 
has a very interesting and a distinct feature where light to moderately heavy (~ 1 TeV) top squarks with negligible 
L-R mixing can be compatible with observed mass of the Higgs boson. Similar possibilities are discussed earlier in 


761 and more recently, in a specific SUSY scenario like the Next-to-Minimal SUSY extension of the 
77 ■ 


generic setups 
SM (NMSSM) 

Taking advantage of this situation, we explore relatively light top squarks satisfying the relevant direct search 
constraints. We further note that a top squark within a narrow mass-window between 197 GeV and 205 GeV f78j , 
which is still allowed by data from the 8 TeV run of the LHC, can be obtained in our scenario for A 5 > 1.2. Such a 
mass-window with > rrit is known in the literature as the ‘stealth’ window [79l - l83| for the top squark when the 
carriers of MET are extremely light (as with the neutrinos in the cascade of the SM top quark). This renders the signal 
characteristics of the top squark rather similar to that of the top quark thus making the presence of the former in 
the data difficult to recognise. Given that the model we consider indeed bears super-light carriers of MET (as we will 
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discuss in section flV All , such a ‘stealth’ scenario can be easily accommodated in our scenarios. Precise measurements 
of the SM top quark properties such as its pair-production cross-section could offer probe to the ‘stealth’ top squark 
scenario. Equivalently, presence of a light top squark may also have some impact on the measurements of the top 
quark mass [8J|. We briefly touch upon the situation with such ‘stealth’ top squarks in section IVIIII 


IV. THE ELECTROWEAK GAUGINOS 

In this section we discuss issues pertaining to neutrino masses and the fermionic sector of the scenario. These will 
help us understand the effect of large neutrino Yukawa coupling (‘/’) in the present context. The electroweak gauginos 
are comprised of the neutralinos and charginos. However, contrary to the MSSM scenario, the neutralinos are Dirac 
fermions in the i?-preserving case. In the i?-breaking case they take a pseudo Dirac form, in general. 


A. The neutralino sector 

The decay branching fractions in various available final states of the top squarks depend crucially on the neutralino 
sector. The neutralino sector in this model differs from that of the MSSM due to the presence of additional fermionic 
fields such as S,T°,R®. In addition, due to the non-vanishing sneutrino vev, the active neutrino mixes with the 
neutralinos. Let us first discuss the neutralino mass matrix in the i?,-preserving scenario, although, ultimately we 
carry out our analysis by considering a mild i?-symmetry violation. This opens up several new and interesting 
phenomenological issues, which we will mention in due course. 


1. The R-conserving case 


The part of the Lagrangian that corresponds to the neutral fermion mass matrix is given by C = (ip 0+ ) T M(P (ip° ) 


where i/j 0+ = (b°,w °, R° d , N c ) and i/j° = (S,T°, H®,v e ). The superscript (±) indicates the respective I?-charges 
which are +1 and — 1. The neutral Dirac fermion mass matrix Mf 1 is given by (52[ 


M x = 


/ M? 

0 

g v u 

V 2 

g Va \ 
V2 

0 

M? 

gv u 

V2 

gva 

V2 

A sv u 

A tv u 

Hu + A svs + A tvt 

0 

V Mr 

0 

-fVa 

fv u / 


( 10 ) 


The Dirac neutralino mass matrix can be diagonalised by a bi-unitary transformation involving two unitary matrices 


V N and U N . The resulting four Dirac mass eigenstates are = 


€ 


with i = 1, 2, 3,4 and Vv = 


Vi = U[j 4’j • With certain simplifying assumptions [52| and with 

A t = As t&n 6w, 

\/2 fMP tan /3 

m r = \ \ — 

g tan 0 w 

the expression for the mass of the lightest neutralino state (x?), i.e., the Dirac neutrino reduces to 


(ID 


= 


v 3 fg sin /3 


V%(Tu + A svs + A tVt) 


A t 


) 

MfMf 


( 12 ) 


Note that only v e acquires a mass at the tree level since only the electron type sneutrino gets a non-zero vev. By 
choosing nearly degenerate Dirac masses for the electroweak gauginos, i.e., ( M— M ~ 0.1 GeV) one can find from 
equation (ED that the Dirac neutrino mass can be in the right ballpark of 0.1 eV even when / ~ 0(1) and assuming 
Mi , M 2 and g. u to be close to a few hundred GeV. The requirement of such a degeneracy between the Dirac masses 
of the electroweak gauginos, however, could be relaxed if one can consider an appropriately small A t ~ 10~ 6 [52j |. 
However, a small A t is not so interesting for our purpose since this results in a diminished contribution to rrih and 
thus, brings back the scenario with multi-TeV top squarks to have the Higgs boson mass at the right ballpark. At 
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colliders, a direct attempt to probe this connection would inevitably involve the heavy right handed neutrino (Nr) 
coupling to an active neutrino ( v e ) and the SM-like Higgs boson. However, while a larger value of ‘/’ enhances this 
coupling, this also pushes up the mass of the right handed neutrino, as conspicuous from eq. m We find that it is 
rather difficult to obtain nriN R < 1.5 TeV in a consistent manner while keeping / ~ 0(1). An immediate probe to 
this at the LHC can be the electroweak process pp —> NRe + h followed by Nr decaying to hv e thus giving rise to a 
final state u e e + hh. Our preliminary study reveals that the corresponding rate could barely reach an attobarn level at 
the 13 TeV run of the LHC and thus one needs to wait for a very high integrated luminosity. A detailed analysis in 
this context is beyond the scope of this paper and we postpone it for a future work (85| . 


2. The R-breaking case 

In the context of a supergravity theory broken spontaneously in the hidden sector, we consider I?-symmetry to 
be broken by a non-zero gravitino mass. The I?-breaking information has to be communicated from the hidden 
sector to the visible sector. We choose anomaly mediated supersymmetry breaking (AMSB) to play the role of a 
messenger |39j, 52]. The I?-breaking Lagrangian contains the following terms 

= Mib°b° + M 2 w°w° -f M$gg + 'y^ A l b L a L b E b + 'y^ A k L a Q k D k 

b- 2,3 fc—1,2,3 

+ E \^ 3k L 2 L 3 E c k + E A bjk L b QjDl + AA H u L a N c + H U QA U U C , 

fc=1,2,3 ?,fc=l,2,3;6=2,3 

(13) 

where, Mi, M 2 and M 3 are the Majorana masses corresponding to U( 1), SU( 2) and SU(3) gauginos, respectively. 
The ‘A’ terms are the standard trilinear scalar couplings. For the rest of this work, we will consider the I?-breaking 
effects to be small, parametrised in terms of the gravitino mass, considered to be roughly around ~ 20 GeV. 

The Majorana neutralino mass matrix containing JS-breaking effects can be written in the basis ip° = 
(b°, S,w°,T, R d: H°, N c , v e ) T as 


qyr = i(^°) T Mf v - 0 + h.c. 

The symmetric 8 x 8 neutralino mass matrix M^ 1 is given by 


/ 

Mi 

Mf 

0 

0 

0 

g v u 

Vi 

0 

g Va 

Vi 

\ 


Mj° 

0 

0 

0 

A sv u 

0 

Mr 

0 



0 

0 

m 2 

M 2 d 

0 

gv u 

Vi 

0 

gv a 

Vi 



0 

0 

m 2 

0 

A tv u 

0 

0 

0 



0 

A sv u 

0 

A tv u 

0 

Tu + A svs + A tvt 

0 

0 



g'v u 

V2 

0 

gvu 

Vi 

0 

Tu + A svs + A tvt 

0 

- fv a 

0 



0 

Mr 

0 

0 

0 

- fv a 

0 

~fv u 


V 

g'v a 

Vi 

0 

gVa 

Vi 

0 

0 

0 

-fv u 

0 

/ 


This can be diagonalised by a unitary transformation 

N*M™ V+ = (M x ) diag . 

We define the two-component mass eigenstates as, 

Xi=Niji>j, i,j = 1, —, 8 . 

Finally the four-component Majorana spinors in terms of the two-component states are defined as 



(14) 


(15) 


(16) 


(17) 


i = 1 , ... 8 . 


(18) 
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The lightest state (xj), which is the neutrino, becomes a Majorana particle. Similar to the Dirac case, the lightest 
eigenvalue of this Majorana neutralino mass matrix resembles the mass of the active neutrino. Other two active 
neutrinos remain massless at this stage. Using the relationships between A s and A t as well as Mr and / as shown in 
section ITV A 11 the mass of the active neutrino can be expressed as [52 ] 


(m^)Tree 


2 [ffA T w 2 (Mf - Mf) sin/3] 2 
V [Mi a 2 + M 2 S 2 ) 


(19) 


where 


a 

5 


2+ A svs + \tVt) tan /3 

g tan 0 W 

V2M®v 2 A t tan /?. 


V2v 2 Xs tan/3(Mj° sin 2 (3 + M 2 cos 2 /?), 


( 20 ) 


It is noteworthy that the parameter A t not only gives a tree level mass to the neutrinos but also help lift the Higgs 
mass through quartic terms generated at one loop level. For example, we have observed that the tree level Majorana 
mass of the active neutrino varies from 0.09 eV to 0.23 eV when A t varies in the range 0.8-1.3, used in Fig. |T| 
which gives correct SM-likc Higgs mass (~ 125 GeV) when top squarks are also not too heavy. Also note that the 
neutrino Majorana mass given in eq. (fl9l) does not depend on the neutrino Yukawa coupling ‘/’. This is because the 
expression has the functional form Mr/ f . Since Mr ~ / therefore, the neutrino Majorana mass becomes devoid of 
‘/’. Similar to the case of Dirac neutrino mass, an appropriately small Majorana mass of the active neutrino requires 
highly degenerate Dirac gaugino masses. Note that (equation ltl9l) j the neutrino Majorana mass is independent of the 
neutrino Yukawa coupling ‘/’. This is in clear contrast to the Dirac case. However, by choosing appropriate values 
of other parameters one can obtain a light active neutrino (x^) with mass 0.1 eV. Nevertheless, some interesting 
observations can be made for various sizes of ‘/’. For example, / ~ 0(1O~ 4 ) gives a sterile neutrino with mass around 
a few keV [58| which can be accommodated as a dark matter candidate. On the other hand, for / ~ 0(1), where we 
obtain a large tree level correction to the Higgs boson mass for low values of tan/3, a light bino-like neutralino (x°) 
with mass around a few hundred MeV (mass of this neutralino is mostly controlled by the //-breaking parameter Mi). 
This MeV neutralino LSP could decay to SM fermions via //-parity violating modes. The probable decay modes could 
be X 2 —> qqv, e + e~v, vvv, qq'e ~, where q, q' are the SM light quark states from the first two generations. However, as 
these involve very small couplings thus resulting in small total decay widths, the decay lengths happen to be much 
larger than the collider dimension T2\. As a result, the LSP neutralino contributes to MET signals. Furthermore, the 
gravitino NLSP (~ 20 GeV) would decay to a photon and the bino-like LSP neutralino. This affects the light element 
abundances which is strongly constrained observationally and results in an upper bound on the reheating temperature 
of the universe Tr < 10 8 GeV [86]. In our scenario, one of the active neutrinos acquire a mass at the tree level. In 
addition, there are one-loop contributions to the neutrino Majorana mass matrix. We observe, that in this model and 
for our benchmark points, the one-loop contributions involving the b — b, t — r loop interfere destructively with the 
Higgs-neutralino loop resulting in a somewhat relaxed bound (~ 20 GeV) on the gravitino mass. This cancellation 
occurs for (m^n where the Higgs-neutralino loop is present. For other elements in the neutrino mass matrix we 
obtain a bound on the relevant RPV operators [52] . This cancellation was not considered previously fU, [HI and as a 
result a stronger bound on gravitino mass was obtained. We have also implemented the model in SARAH-4.4.1, which 
performs a full one-loop correction to the neutralino mass matrix, and we have cross-checked with the spectrum file 
that gravitino mass in the ballpark of 0(10) GeV is consistent with light neutrino masses and mixing. 


B. The chargino sector 


Just like the neutrino and the neutralinos would mix, the charged lepton mixes 3 with the charginos. This results 
in an extended chargino mass matrix [63] compared to the MSSM. The relevant Lagrangian after //-breaking consists 


3 The mixing between the charged leptons and the charginos gives rise to deviation of the Z to charged lepton couplings. Such a deviation 
is very much constrained from the electroweak precision measurements leading to a lower bound on tan/3 >2.7 [53 . l45j| . It is also 
pertinent to mention that an upper bound on tan/3 comes from r Yukawa coupling contributing to the ratio R T = T(r —»• ev e i> T ) /T(r —>• 
[3l l45ll . Choosing to be around 1 TeV corresponds to tan/3 < 70. 
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of the following terms 

C ch = M 2 w + w~ + (Mf - gv T ) T~w + + V2 X t v u T+R~ + gv u H+ w~ - g u H+ 

+ A t v t H+R-[ - A s v s H+Rj + gv a w + e~[ + ( M 2 + gv T ) T+w~ + m e e R e~ L + h.c. 

(21) 

The 4x4 chargino mass matrix written in the basis ip Y = (w + , T+, H+, e R ) T and ?/>“ = (w~ , T)~, R^ , e£) T is given 
by .' 

( M 2 Mf - gv T 0 gv a \ 

M + gvr 0 \/2v u X t 0 ^ 2 ) 

gv u 0 -g u - X s v s + X T v T 0 

\ 0 0 0 TTl e ) 

Again, this can be diagonalised by a bi-unitary transformation UM C V T = M|j. The chargino mass eigenstates 
(two-component) are written in terms of the gauge eigenstates in a compact form as 

XT = UijipJ, 

xt = Vi^j. (23) 

The four-component Dirac spinors written in terms of the two-component spinors take the form 

#=^J> (i = 1,-., 4). (24) 


The lightest chargino (Xi) corresponds to the electron, x 2 is the lightest chargino state reminiscent of the lighter 
chargino in the MSSM with mass of 0(100) GeV. It is also pertinent to mention that i/)~ and ipf would also include 
g~[, t l an b M/ji t ri respectively. However, as discussed in section [TT] only the electron type sneutrino acquires a vev, 
and the vev of the other two sneutrinos can be rotated away without any loss of generality. Therefore, g and r do 
not mix with the chargino states. 


V. THE TOP SQUARK SECTOR 


In this work we concentrate on the third generation squarks, mainly the top squarks, which play important roles in 
lifting the Higgs boson mass. Scenarios with light top squarks draw their motivations from the ‘naturalness’ argument. 
They also provide rich and interesting collider signatures. As discussed in section Mil the model which we consider 
here gives us the opportunity to study such light top squarks. Furthermore, l?-symmetry prohibits any trilinear scalar 
couplings (the ‘A’ terms) and Higgsino mass parameter (the g term). Therefore, we investigate a situation where 
both the top squarks are light (~500 GeV) and have negligible chiral mixing, which originates from small f?-breaking. 

The relevant terms in the top squark mass matrix are generated from the f-term, the D-term and the soft terms. 
The SU(2)l and XJ{ l)y contributions to the JD-fields are given by 


D a 

D y 


9 


H\j a H u + L\r a U + q\ L T a q iL + T^tI + y/2 \M°T a + Mf T at 


1 


:9 


2 

y/2M? 


q i L 

~+~ 1 j. _ 4 . 2 -a. ~ ' 

H' U H U - L\Li + 2e* iR e iR + -ql L q iL ~ ^u} R u iR + -d\ R d iR 


S + S f 


(25) 


The r and A matrices are the generators of SU( 2)^ group in the fundamental and adjoint representation, respectively. 
From equation (1251) it is straightforward to calculate the elements of the mass-squared matrix in the top squark sector, 
which in the basis (t R ,t R ) turn out to be 


{M~) n = m | 3 + to 2 + m| cos 2 $ Q - | sin2 ^ v s M f + y/2gM 2 v T , 

(M?) 12 = (M?) 21 = 0, 

2 4\/2 

(M ~) 2 2 = ml 3 + m 2 + -to| sin 2 6 W cos2/3- —g'M°v S - 


( 26 ) 
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Note that in absence of ‘ A ’ terms and the p term, the off-diagonal entries vanish and hence the top squark sector 
is devoid of any chiral mixing. Thus, the left- and right-chiral states are equivalent to the mass eigenstates. Such 
a ‘zero’ mixing situation can be contrasted with the MSSM, in which, a substantial mixing is generally required to 
obtain the observed value of the Higgs mass. Motivated by the recent collider bounds on the masses of the top squarks 
(depending on its various decay modes), to be justified in some detail in section IVCl we choose t\ ~ and t 2 ~ tz,. 

At the LHC, top squarks are being searched in their direct production, pp —► tt*, followed by their subsequent 
decays in various possible modes. Out of these, the decays that are relevant to our scenario |87H90| are the R-parity 
conserving ones 


t -A bx + and (27) 

and the modes that violate R-parity when a top squark could decay to a bottom quark and a charged lepton [in], |)2j 

t^M + . (28) 

These channels are of major relevance in the context of our model. The reasons are twofold: first, in the large 
‘/’ scenario, we obtain a light bino-like neutralino with mass around a few hundred MeV, in addition to an active 
neutrino. Therefore, the channels with top squark decaying to a top quark and a bino-like neutralino and/or an active 
neutrino open up. Secondly, top squark decaying to a bottom quark and a chargino is also important. Additionally, 
top squark decaying to a bottom quark and an electron becomes an interesting channel to look for. This decay mode 
is predominantly controlled by the R-parity violating operator A' 133 . In the framework of the MSSM with R-parity 
violation, strong limit on this particular coupling exists from the neutrino Majorana mass, |A' 133 | ~ 3.4 x 10” 3 

Hence, in such a scenario, the resulting decay rate becomes highly suppressed. However, in the present context, A' 133 
is identified with the bottom Yukawa coupling yb . The smallness of the neutrino mass is then explained through small 
R-breaking effect, parametrised in terms of a small gravitino mass [52| . Thus, a large decay rate for t —> be + becomes 
a generic feature in our model. In principle, t could also decay to bp + and br + via RPV couplings A 2 33 and A 333 . 
However, these decays are subdominant compared to t —> be + because of the stringent constraints on the relevant 
couplings as discussed later. We note in passing that t —> tG is also a possibility but highly suppressed 0 for a 
gravitino of mass ~ 20 GeV in the present context. 

The relevant Lagrangians are worked out in the four component notation following [95j, j9f| and are given by 



C 


o 


= —t 


h.c 


ytPLNiQ + —^=| gPRN a + y PrNu 


tLXi +1 


3\/2 


PlNh - ytPnNiQ 


tRXi 


(29) 


and 




\ 

1 


C tb^ = b 

a? 

1 

tLXi + b 

CO 


tRXi + A/ 33 tLPLUi4Xib + h.c., 


(30) 


where, A' 133 = yb = yyjjbj; the bottom Yukawa coupling and yt = the top Yukawa coupling and m* and mb 

are the top and the bottom quark masses, respectively. The neutralino and the chargino mixing matrices N. t j. Uij 
and Vij are as defined earlier. Note that, for i = 1 (corresponding to \i — e ~) the mixing matrix elements U\\ and 
V 13 are suppressed. In the following subsections we briefly discuss the salient decay modes of the lighter (t± ~ £#) 
and the heavier ( t 2 ~ £l) top squarks. 


A. Decay rates of ti (as to) 
The partial decay widths of t\ in the ty 1 / and b\f inodes are given by 

1 


T(fi -A t\i ) — 


167TTO- 


V Ri + C Ri j ( ~ m t - ) - ^VRiCmmtm^ 


4 i 4 i 4 ci 2 Z n 2 2 ci ‘2 2 

m~ + m~o + m f — 2m~ m f — 2m~om f — 2m~ 0 m 7 

^ 1 11 L Xi 1 


Xi 


Xi 1 1 


(31) 
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and 


where, 


r (*i bxf) = 


167r?n~ 


l Ri 


Pm 


2 2 

— m h — ra~+ 
0 xT 


- 4:amPmmbm~+ 


a 

■ m~+ 
x7 


mi — 2 m~ ml — 2m^+ml — 2mi+m~ 
0 1 1 0 xj 0 xj ( i 


,s ‘ = 

C Hi = 

OiRi = 0 , 

P Ri Vt^i3- 


(32) 


(33) 


We note down a few important observations below. 

• In the large ‘/’ case we obtain a light (~ few hundred MeV, governed by the .R-breaking Majorana mass Mi) 
bino-like neutralino (y!))- This is because of the presence of MrN c S term in the Lagrangian, where the coefficient 
Mr becomes very large (~ 10 5 GeV) for an order one ‘/’. This results in forming a heavy pseudo-Dirac pair 
with mass ~ Mr and makes the lightest eigenvalue very small and predominantly bino-like. 

• The Dirac wino mass M® is considered to be heavy to evade bounds from Z boson coupling to electrons fibl ] 
(see also section HV Bl) . The [i u parameter, which controls the mass of the Higgsino (both neutral and charged) 
can vary between the electroweak scale (~ 200 GeV) and a much larger value, i.e., a few TeV. 

• Based on the above discussion and with the help of equations (l29l) - (l33ll . we find that t\ would decay into t \° 3 4 

and b \2 ■ The neutralino can be both bino- or Higgsino-like whereas the chargino would only be Higgsino-like, 
assuming the Higgsino mass parameter rn ^) << M®. 

• We expect the dominant decay modes of t\ to have the Higgsino-like neutralinos or chargino as the decay 
products rather than the bino-like neutralino. This is because of the enlarged couplings for the former which 
are proportional to the top Yukawa coupling, y t . 

• Finally, in the limit when , the top squark cannot decay to an on-shell top quark and a Higgsino- 

like neutralino or a bottom quark and a Higgsino-like chargino due to phase space constraints. Therefore, the 
dominant channel would only be t\ —> fy®, where y!j is the bino-like MeV neutralino. Moreover, t\ —> tv e would 
also contribute to MET, although the branching is suppressed due to the small neutralino-neutrino mixing. 


B. Decay rates of £2 (~ tr.) 
The partial decay widths of £2 are given by 

1 


and 


T (£2 — >tXi) — 


167rm2 

t 2 


VLi + Cm ) ( m t 2 - m t ~ ) - 4 VLiCurntm. 


4 1 4 1 4 r>22 r> 2 2 q2 2 

m~ + m~o + m f — 2m~ m t — 2m~ 0 m t — 2m~ 0 m~ 

b 12 L Xi L Xi t-2 


t 2 1 Xi 


T (£ 2 ->• bxt) = 


167rm2 
12 . 


a Li + PLi )\ m t 2 ~ m b- ) - 4 a Li p L im b m~+ 


2 
L 

m~ + ml+ + mi — 2m~ ml — 2ml. + ml — 2mt+ m~ 
1 2 x7 0 *2 0 x7 0 x7 *2 


(34) 


( 35 ) 
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where, 


VLi 
C Li 


&Li 


yt N l6 , 

7 ! (^3 + 1^1 
-gUn , 


fai = o. 


(36) 


In addition, for f 2 we also have the interesting possibility of f 2 ~t be + . The corresponding partial decay width is given 
by 


T(f 2 —>• be + ) 


vl l^4| 2 

167T 


"t 2 - 


(37) 


Some features of f 2 decays are as follows: 

• The decay t 2 -» be + is an interesting possibility. This faces no suppression from the phase space and the decay 
rate is proportional to the bottom Yukawa coupling, yb which grows with tan /?. Hence a substantial branching 
fraction in this mode is expected at large tan/3 and for a fixed top squark mass. 

• When fj, u < mj 2 , t 2 would decay to Higgsino-like chargino and neutralinos. Also, decay to a bino-like neutralino 
is a possibility. However, a quick look at the couplings in equation (ERJll would suggest that the decay to Higgsino- 
like neutralinos (x®, xi) is rjLi ~ yt-enhanced and hence, is more probable than a decay to a bino-like neutralino 
(suppressed by g'/ 3 in the coupling) or to a Higgsino-like chargino (suppressed by g times the wino component 
of the lighter chargino, Un). 

• Again, for g u > mj 2 , decays of t 2 to Higgsino-like neutralinos and charginos are kinematically barred. Under such 
a circumstance, t - 2 mainly decays to a bottom quark and an electron (positron). The decay mode t 2 (tL) —> t \ 2 
is again suppressed because of a (comparatively) small involved coupling. Finally, am = /3m = 0 reflects the 
absence of Hd in the Lagrangian, which has been integrated out from the theory. 


C. Bounds on top squarks 


• Recently ATLAS measured the spin correlation in the top-antitop quark events and searched for top squark pair 
production 97| in the pp collisions at -Js = 8 TeV centre of mass energy and integrated luminosity (£) of 20.3 
fb -1 . This particular search has ruled out top squarks with masses between the top quark mass and 191 GeV 


with 95% confidence level. A very recent study [78| reveals that the window of 197 GeV < mq < 205 GeV, in 
the so-called ‘stealth’ regime (i.e., with vanishing LSP mass), cannot yet be ruled out. 


• Dedicated searches for pair-produced top squarks decaying 100% of the time to bottom quarks and lighter 
charginos have been performed [98;, ;99| within the framework of the MSSM. For a chargino with mass close to 
200 GeV, the top squark below 470 GeV has been ruled out at the 7s = 8 TeV run of the LHC. In our scenario, 
tn decays to this particular channel if g u < nij R . Although in our model the corresponding branching fraction 
is less than 100%, we take a conservative approach and respect this bound. In addition, this search gives the 
most relaxed bound on the mass of the top squark which is relevant to our analysis. Hence in the present study 
we choose to be the lighter top squark, i.e., t\ ss t r. 

• Another decay mode of the top squark relevant to our scenario is t —> t\i, where implies the lightest 
supersymmetric particle in the MSSM (in our work, however, Tj 1 is identified with the active neutrino and \2 
represents the lightest bino-like neutralino). At the 8 TeV run of the LHC, top squark with mass below 550 
GeV is ruled out at 95% confidence level [98, [99j with the assumption BR ( t\ —> fy?) = 100%. This bound 
applies for a massless neutralino (m^o = 0). For heavier neutralinos in the final states, this lower bound on the 
top squark mass can be relaxed further. Note that in the large ‘/’ scenario we find a light super-light bino-like 
neutralino with mass around a few hundred MeV which thus attracts this bound on the mass of the lighter top 
squark. 


• A top squark decaying via i?-parity violating mode has also been probed by the LHC experiments. If a top 
squark undergoes an R -parity violating decay only to a bottom quark and an electron, a stringent lower bound 
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[El [92i l exists on the top squark mass with mj > 900 GeV. 4 Accommodating an even lighter top squark, which 
is central to our present work, thus requires a situation where such a bound is preferentially applicable to the 
heavier top squark state (t 2 ) of the scenario. As described in section IVBl only tL could decay to a bottom quark 
and an electron (positron). Hence we choose the heavier top squark t 2 to be the t^, i.e., t 2 ~ Note again 
that our consideration is pretty conservative and, as we would find in section eh for generic scenarios where 
such a decay canjiave a branching fraction below 50%, the bounds can get considerably weaker thus allowing for 
an even lighter t 2 . Phenomenological discussions on top squarks undergoing such an A.-parity violating decay 
can be found in references UMozl. 


The squarks from the third generation have understandably attracted a lot of attention in the recent times. The 
flavor changing decay of the top squark, t\ —> cx± has been analysed in great detail in :77l [MnH . Recent searches 
performed by both ATLAS and CMS collaborations have looked into this channel extensively and ruled out top squark 
masses below 300 GeV mmm- Top squarks decaying to a top quark along with a neutralino (LSP or NLSP) 
has also been probed in various SUSY models. A lower limit close to 1 TeV for the top squark mass can be obtained 
at ^fs =14 TeV and with the high luminosity option Ill4l4llfi |. In addition, thorough phenomenological studies have 
also been performed in the decay of top squark into a bottom quark and a chargino 117. Ill8j . 

We note in passing that interesting final state signatures can be obtained for the decays of bottom squarks as 
well PI. For example decays to a bottom quark and the bino-like neutralino with a branching ratio close to 71% 


and bji decays to a top quark and an electron with a branching close to 47% in the BP-I scenario. Such a branching 
would imply 2 6-jets+^r or 2 &-jets+4 leptons+,i?T in the final states, respectively. The most stringent limit on the 
mass of the bottom squark comes from the search where it decays to a bottom quark and the lightest neutralino 


(LSP) with BR (b —> &\°)=100%. Bottom squark mass up to 700 GeV has been excluded at 95% confidence level for 
neutralino mass less than 50 GeV Il20| . We note that all the benchmark points are chosen in a way such that they 
satisfy existing bounds on top squark mass under various circumstances pertaining to its decay. 


VI. THE BENCHMARKS AND THE FINAL STATES 


In this section we discuss a few benchmark scenarios that would be broadly representative of the phenomenology 
that is expected of the framework under consideration. As mentioned earlier in section EH we embed the model in 
SARAH (v4.4.1) [7314751 . We use the low energy output of SARAH (v4.4.1) and generate the SUSY spectrum using 


SPheno (v3.3.3) 


1211 122 


energy 

FlavorKit 123] is used to ensure benchmark points are consistent with all rel evan t 
Higgs boson cross-sections and signal strengths are computed using HiggsBounds [l24l — 
As discussed earlier, we will mainly consider two regimes, viz., n u > my. „ and 


flav or violating constraints;_ 

Il27| and HiggsSignals jl28l. 12j.j . J uiovuoocu i.ai , vv... win mainly y.'.nioiuy,i iwu 1 ii, vi/,., / lll, t\ 2 

Hu < . For each case, we point out the dominant decay modes of both t\ and t 2 . These dictate the types of 

interesting signatures at the LHC for each of these cases. 


A. Case 1: //.„ < mq 

Two benchmark points for this case are shown in table El A relatively low value of ji u ( = 200 GeV) is chosen for 
both the benchmark points. The masses of the Higgsino-like chargino and the neutralinos are mainly controlled by 
Hu. We assume the singlet and the triplet vevs to be small; roughly to be around 10 -4 GeV. The Dirac gluino mass 
(MjP) is considered to be 1.5 TeV. Since we are considering a small amount of -R-breaking, the Majorana gaugino 
masses are roughly around a few hundred MeV. Fixing the order parameter of -R-breaking, i.e., the gravitino mass 
0(10 GeV), fixes these soft SUSY breaking parameters. Both As and A t are considered to be large, which for 
large tan/3 (> 20), provide significant radiative corrections to the Higgs boson mass through one loop quartic terms. 
Such a choice allows us to have a situation where both the top squarks are moderately light. Note that the chosen 
values of m^(~ 470GeV) for BP-1 and BP-2 are expected to be consistent with the latest LHC bounds discussed in 


4 Note however, that if t —> br + opens up, the corresponding lower bound on the top squark mass can be relaxed. The decay t —r br + 
is mostly controlled by the R-parity violating coupling A 333 . The existing bound on this particular coupling is much relaxed: Agg 3 < 
1.4 cos/3 [45l | and can be saturated for small values of tan 3 (< 5). However, in the present scenario we confine ourselves in the limit 
where tan/3 20, which renders the decay t — » br + insignificant. On the other hand, the decay t —> hfC is also negligible because of 
the strong constraint | A 233 1 < 6.8 X 10 “ 3 cos/3 llOd . 
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section lYCl The bound assuming BR(fi — > b\t) = 100% is evidently satisfied while the one (mq > 550 GeV) that 
assumes BR(H —> txi) = 100% is not applicable here (see table ITTT1) . Furthermore, the neutrino Yukawa coupling ‘/’ 
is chosen to be 1. Hence, to have the active neutrino mass in the right ballpark, we need to consider the Dirac bino 
and wino masses to be almost degenerate. As conspicuous from equation (flTH) . this degeneracy provides a suitable 
suppression to the neutrino Majorana mass when the Dirac gaugino masses themselves are roughly around a TeV or 
so. Some of the low energy flavor violating branching ratios (which satisfy the respective experimental constraints) 
are also shown in table El 


Parameters 

BP-1 

BP-2 

Mf 

1200 GeV 

800 GeV 

M 2 d 

1200.1 GeV 

800.1 GeV 

Mf 

1500 GeV 

1500 GeV 

f~tu 

200 GeV 

200 GeV 

m 3/2 

20 GeV 

20 GeV 

tan /3 

23 

35 

(urS)33 

2.3x10 s GeV 2 

2.5 x 10 s GeV 2 

i m o) 33 

5.5x10 s GeV 2 

6.1 x 10 s GeV 2 

/ 

1 

1 

Vs 

2 x 10~ 4 GeV 

1.5 x 10 -4 GeV 

Vt 

> 

CD 

O 

1 

O 

10“ 4 GeV 

A s 

1.130 

1.116 

Bfl L 

-(200 GeV) 2 

-(200 GeV) 2 

ts 

(174 GeV) 3 

(174 GeV) 3 

Observables 

BP-1 

BP-2 

m h 

124.9 GeV 

125.7 GeV 

m {l 

566.2 GeV 

580.5 GeV 

m h 

918.5 GeV 

904.8 GeV 

m~o = m Ue 

0.01 eV 

0.04 eV 

m^o (bino-like) 

167.9 MeV 

168.3 MeV 

m-o 

211.5 GeV 

213.8 GeV 


211.5 GeV 

213.8 GeV 

m~+ = m e 

0.51 MeV 

0.51 MeV 

m~+ 

243.8 GeV 

247.1 GeV 

Flavor Observables 

BP-1 

BP-2 

BR (B -> Xsy) 

3.4 x 10 -4 

3.3 x 10" 4 

BR (B° s -7 MM) 

2.4 x 10" 9 

2.5 x 10" 9 

BR(m —> ey) 

3.8 x 10" 24 

4.9 x 10 -24 

BR(/r —» 3e) 

3.0 x 10~ 26 

4.0 x 10 -26 

m 77 

1.05 

1.06 


TABLE II: Benchmark sets of input parameters in the large neutrino Yukawa coupling (‘/’) scenario and the 
resulting mass-values for some relevant excitations for /i u < mq (case 1). M 3 denotes the Dirac gluino mass. Also 
indicated are some of the relevant flavor observables and their values, all of which are currently allowed by 
experiments. The corresponding values of /x 77 (the estimated Higgs di-photon rate compared to its SM expectation) 

are also mentioned. 


• Decay branching fractions of t\ 

The dominant decay branching fractions of t\ for BP-1 and BP-2 are indicated in tabic 1IIII \2 is the Higgsino-like 
chargino and 4 are the Higgsino-like neutralinos. X 2 is the bino-like neutralino with mass in the ballpark of a few 
hundred MeV. As can be seen from table Mil the top squarks, once produced in pairs, can undergo both symmetric as 
well as asymmetric decays. Table HV1 lists all possible final state topologies. However, in the context of this work, we 
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Decay modes 

BR for BP-1 

BR for BP-2 

ti-tbxi 

52.5% 

51.7% 

ti -»• t X3 

20.0% 

20.1% 

ti -¥ t x° 

20.0% 

20.1% 

tl H► t X2 

6.0% 

6.6% 


TABLE III: Decay branching fractions of t\ in BP-1 and BP-2 for p u < 


will mainly consider the dilepton final states accompanied by 6 -jets and MET. Such a final state might arise when the 
top squarks, on being pair-produced, undergo the decay t\ —> bxt ■ xt i n turn, decays to a W + and a x 2 (x?) with a 
branching ratio close to 90% (10%) followed by W -s decaying leptonically. Although a semileptonic {ivjj) final state 
from TE-decays is a good compromise between the rate and the cleanliness of the signal, we go for a cleaner channel 
where both the W bosons decay leptonically. As shown in figure EE such a topology leads to a final state 2 6 -jets + 2 
leptons +J&f. A similar final state could also arise from other decays of the lighter top squark, such as those involving 
f i —f tx 2 as shown in figures [2a] and [2b] However, the effective branching ratio is rather suppressed. 

We note in passing that the various decay combinations shown in table Mil could also provide exotic multilepton 
and multijet final states depending on the leptonic or hadronic decays of both W± or Z boson. For example, decays 
such as t\ —> txy 4 could give rise to a 2 6 -jets + 6 leptons + final state. Some relevant final states arising from 

the decays of t\ are tabulated in table IIVI The branching fractions in BP-1 and BP-2 are rather similar since we 
are dealing with similar top squark masses. Also, top Yukawa coupling is practically insensitive to larger values of 
tan/3, as considered in our study. This results in similar branching fractions in the bx + mode in BP-1 and BP-2. The 
dynamics of other decays are essentially controlled by the gauge couplings and therefore, they remain similar. 


p u < : Decays of f 1 

pp-> titj bW^cl —>■ 2b + 2W+M 

pp tifM- tX.2 ix°^bW+5° bW^g ->• 2b + 2W+Mf 
pp ->■ Fit* -> bxj + h.c. -> bW+^l bW^o + h.c. -* 26 + 2W + 

pp -¥ titl tx 3/4 tx 2 + h.c. —► bW+Zx .2 bW~x .2 + h.c. —¥ 26 + 2 W + Z + JiY 
pp ->■ ht* -> b& t^ i + h.c. -» bW+^l bW^xl + h.c. -> 26 + 2W + Z + 
pp -» Ut\ ->■ tx ° 3/4 ixl / 4 -» bW+Zx° 2 bW~Zx° 2 ^2b + 2W + 2Z 

TABLE IV: Possible final states arising out of various decay modes of t\ when fi u < rri^. 


• Decay branching fractions of t 2 


The dominant decay branching fractions of t 2 for BP-1 and BP-2 are shown in table [Vj The pattern can be justified 
following the discussion in section fV Bl The branching fractions of t 2 to the three modes indicated are comparable. 


Decay modes 

BR for BP-1 

BR for BP-2 

£2 —^ 6 

27.8% 

47.2% 

t 2 ~¥ t X3 

35.7% 

26.0% 

*2 —> t X4 

35.7% 

26.0% 


TABLE V: Decay branching fractions of t 2 in BP-1 and BP-2 for which p u < mj 2 . 

These lead to distinct final state signatures with appreciable strength. Possible final states arising from the decays 
of t 2 are listed in table ED The decay channels t 2 —» tiZ(h) are absent due to negligibly small mixing between the 
left and the right chiral states of the top squark. A remarkable point to note here is the significant decay branching 
fraction of t 2 to a bottom quark and an electron. To reiterate, this decay rate is proportional to A' 133) which is 
identified with y^. Therefore, the corresponding decay rate is large. Also, because of the large difference between 
the mass of the decaying particle ( 62 ) and the total mass of the particles in the final state (mb + m e ), the final state 
electron is expected to be hard. The schematic diagram for such a process is presented in figure [3] Other decay 
modes, presented in table ED are similar to the previous case where the decay products of t 2 are a top quark and a 
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(a) Asymmetric decays of the top squarks with one decaying 
to bxt while the other decaying to thus leading to 
2 b + 2 £ +JZt final state. 


(b) Symmetric decays of top squarks with both decaying to 
t\ 2 , thus leading to 2 b + 21 +J£t final state. 




(c) Symmetric decays of top squarks with both decaying to 
bx ^, thus leading to 2 b + 2 £ +JZt final state. 

FIG. 2: Final state topologies yielding 2b + 2£ + in the 


fiu < Tn ^2 : Decays of t 2 
pp —1 12^2 ~t be+ be~ —>■ 2 b + e + e - 

pp —► t 2<2 —> be+ t%° 3 / 4 + h.c. — > be+ bW~Z\ 2 + h.c. —> 2 b + W + Z + e + + JSY 
PP ^ t^Q^bW+Zx° bW~ Zx° —> 26 + 2W + 2Z + J&r 


TABLE VI: Possible final states arising out of various decay inodes of 1 2 when p, u < mt 2 ■ 


Higgsino-like neutralino. These would further decay to give a final state comprising of 6 leptons+2 6 -jets +Xt- The 
variations in the branching fractions as we go from BP-1 to BP-2, as can be seen in table El are due to changing 
bottom Yukawa coupling as tan /3 changes. 


B. Case 2: p u > 2 

As opposed to the previous case, we consider the situation where n u > a . To have the Higgs boson mass in 
the right range we tweak A s- The soft masses (m u )| 3 and (toq) g 3 are modified to get different top squark masses 
satisfying relevant LHC constraints, vt just takes a different in sign when compared to BP-1 and BP-2 only to exclude 
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FIG. 3: 2 b+e + e final state arising from both f 2 -s decaying directly to a bottom (anti-bottom) quark and a positron 

(electron). 


tachyonic states. All the other parameters are kept fixed to their values in table ITT] Due to such a choice of /j, u , a 
top squark cannot decay to Higgsino-like chargino and neutralinos. Note that, in the present case as well, the masses 
of the top squarks are chosen in such a manner that they satisfy the present experimental bounds. Values of flavor 
observables are checked to satisfy experimental constraints. However, those are not shown explicitly this time. 


Parameters 

BP-3 

BP-4 

fJ-u 

1500 GeV 

1100 GeV 

vt 

-10“ 4 GeV 

-10“ 4 GeV 

(™«) 33 

4xl0 5 (GeV) 2 

5x10 s (GeV) 2 

(™q)s 3 

5.2x10 s (GeV) 2 

6x10 s (GeV) 2 

As 

1.09 

1.06 

Observables 

BP-3 

BP-4 

m h 

126.6 GeV 

126.1 GeV 

m h 

728.7 GeV 

802.8 GeV 

m i 2 

909.5 GeV 

908.8 GeV 

m~o = m Ue 

0.03 eV 

0.12 eV 

m^o (bino-like) 

175.8 MeV 

175.6 MeV 

m~o 

1202.1 GeV 

804.3 GeV 

m xZ 

1202.2 GeV 

804.3 GeV 

m~+ = m e 

0.51 MeV 

0.51 MeV 

m~+ 

_xl_ 

1304.2 GeV 

877.4 GeV 


TABLE VII: Same as in U but for an extra sign on vt and for BP-3 and BP-4 for both of which /i u > mj (case 
2). Values of flavor observables, not shown here explicitly, satisfy all the experimental constraints. 


• Decay branching fractions of t\ 

The decay branching fractions for t\ when fx u > 2 are shown in tabic IVIIII for the benchmark points BP-3 and 

BP-4. An interesting point to note here is that the top squark decays only to a top quark accompanied either by a X 2 
( figure l2bl) or a v e (figure 0} both of which are carriers of MET while the former being the dominant one. Again, both 


Decay modes 

BR for BP-3 

BR for BP-4 

ti tx 2 

ti —> tv e 

87.8% 

12.2% 

94.6% 

5.3% 


TABLE VIII: Decay branching fractions of t\ in BP-3 and BP-4 for which /i u > . 

symmetric and asymmetric decays of the pair produced t±-s are possible. These would lead to 2 b- jets+2£ + $t final 
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states. Note that more exotic final states with a larger lepton multiplicity would be absent as heavier Higgsino-like 
neutralino(s) will now be missing in the cascades of t\. This is in sharp contrast with what is expected for p, u < mj 
(case 1) as discussed in section lYI Al and thus, may be exploited to distinguish between these two broad scenarios. 

• Decay branching fractions of £2 

Similarly, the absence of a light Higgsino-like chargino and neutralinos implies £2 would dominantly decay to a 
bottom quark and an electron (positron). The branching fractions of £2 under such a circumstance are presented in 
table HXl 


Decay modes 

BR for BP-3 

BR for BP-4 

i 2 —> b e + 
t 2 “S' t X2 

97.8% 

1.8% 

98.7% 

0.9% 


TABLE IX: Decay branching fractions of £2 for BP-3 and BP-4 for which p u > mj 


VII. COLLIDER (LHC) ANALYSIS 

In this section we present the setup and the results of the simulation we carry out at the 13 TeV LHC for pair- 
produced top squarks that eventually cascade to the final states discussed in section IVT1 


A. The simulation setup and reconstructing the physics objects 


We have implemented the model in MadGraph5_aMC@NL0 [1301 ]. Events for both signals and backgrounds are gener¬ 
ated using the same. We use the parton distribution function CTEQ6L1 |l3l| and a factorisation/renormalisation scale 
of , yrnf^rn^ for generating events at the lowest order. The inclusive rates are then normalised to their respective 
values obtained after higher order corrections as given by MadGraph5_aMC@NL0 in the cases for the SM background 
and Prospino2 (v2.1) damn for the case of the SUSY productions. Appropriate branching fractions are ob¬ 
tained from the spectrum generator SPheno 121] which, in the first place, is generated by SARAH |73l - l75| . We note in 
passing that the production cross-section for the top squarks (at the tree level) in this model is same as in the MSSM, 
considering only the dominant strong interaction. 

Events in the LHE format are fed into Pythia-6 .4.28 Il34[ for showering, hadronisation and jet formation. Clus¬ 
tering of jets is performed with the built-in Pythia module PYCELL which employs a cone algorithm and incorporates 
appropriate smearing of the momenta. In PYCELL we allowed for an angular coverage of |?;| < 5 for the hadron calorime¬ 
ter with a cell-segmentation of Ap x A</> = 0.1 x 0.1 which resembles a generic LHC detector. A cell is required to have 
a minimum value of deposited Et = 1 GeV for it to be considered. A jet-cone radius of A R(j,j) = 0.4 is employed 
for finding jets. A minimum summed Et of 20 GeV is required within such a geometry for the configuration to be 
considered as a jet. Ultimately, formed jets within |? 7 | < 2.5 are considered in our analysis. Care has been taken to 
isolate final state leptons by imposing the following cuts and isolation criteria: 


• To select leptonic events we have used pj, > 10 GeV and \pe\ < 2.4. 

• Lepton-lepton separation has been done by choosing A R(£',£) > 0.2, where A R = (A 77) 2 + (A (j>) 2 . 

• Subsequently, to separate leptons from jets we have used A R(j, t) > 0.5. 

• Finally, the sum of the energy deposits of the hadrons which fall within a cone of A R < 0.2 around a lepton, 
must be less than 10 GeV. 


In this work, by leptons we mean only electrons and muons for which the detection efficiencies are generally very 
high, unlike the r lepton. We have used a minimum pr cut of 10 GeV and 17 GeV to isolate muons and electrons, 
respectively. To estimate the number of 6 -jets in the final state, a flat (but somewhat conservative) 6 -tagging efficiency 
of 60% has been used. 
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B. Top squark pair-production cross section 


The phenomenology we discuss in this work crucially depends on the rate of pair-production of the top squarks. It 
is to be noted that at the lowest order these rates are the same as in the MSSM. Considering the dominant strong 
contributions in the rates, the variation of the same is only dependent on the mass of the top squark, irrespective of 
its chiral content. As a quick reckoner for this basic rate, we present the same as a function of mz in figure [5] for the 
13 TeV run of the LHC. Appropriate A'-factors as obtained from the package Prospino2 |l32l . ll33f are already folded 
in. To this end, we fix all the parameters as given in BP-4 except for the right handed soft squark mass, which we 
varied from —4 x 10 4 GeV 2 < ( 771^)33 < 8 x 10 5 GeV 2 . Such a choice would surely move the Higgs mass away from 
the allowed range. However, we are here merely concentrating on study of the production cross-section for tt*. The 
parameters such as A 5 , A t can be adjusted to fit the Higgs mass, which is unlikely to affect the production rate any 
significantly (via unknown higher order effects). 




mj{ GeV) 


FIG. 5: Variation of production cross section for a pair of top squarks at the 13 TeV LHC. Other parameters are 
kept fixed at values given for BP-4 in tabic I V I FI (see text). Rates include appropriate A'-factors obtained from the 

package Prospino2. 
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C. The Standard Model backgrounds 

As mentioned earlier, we would mostly concentrate on the final states with 2 6 -jets +2 leptons+Xr and 2 b- 
jets+e + e _ . In the first case, we consider only the most dominant background coming from tt production when 
both the top quarks decay leptonically. In order to have a realistic normalisation of this background, tt events gen¬ 
erated at the lowest order (LO) using MadGraph5_aMC@NL0 and the size of the event sample is scaled appropriately to 
correspond to its next-to-leading order (NLO) + next-to-next-to-leading log (NNLL) cross-section (« 816 pb) [l35| . 
The heavier top squark undergoes a significant decay to a bottom quark and an electron. The dominant background 
comes from the direct production of a pair of bottom quarks with one of them radiating a Z or 7 * which subsequently 
produces a pair of e + e _ . This background can be largely suppressed by putting an on-shell Z veto for the e + e - 
pair. To be conservative and for the robustness of the estimate, the NLO computation is done with two 

additional jets (10 GeV < < 60 GeV) in the final state. The SM cross-section for pp —> 6 &e + e _ +jets we used is 

9.43 pb. 


D. Event selection 

To optimise the signal to background ratios, we now have to adopt a set of event selection criteria. Towards this, 
various appropriate kinematic distributions are studied for both signals and the backgrounds. We present our study for 
two broad scenarios discussed in section ED i.e, for p u < nip and p, u > rrip . For each of these cases, two different 

final states are considered, viz., 2 &-jets +2 leptons +Xt and bbe + e~, arising from t± and t 2 decays, respectively. 


1. Case 1: fi u < nip 2 


Here we discuss the decays of both the top squarks pertaining to the case where p u < nip . 

• pp —> tit* —>• b\ 2 /t \§ —¥ 2 &-jets + 2 leptons + Xr (figure [ 2 ]) 

Such a final state could arise from top squarks decaying to b\f and/or tx 2 - The final state leptons arise from 
the decays of W bosons. In addition, xt is somewhat heavier than the top quark for both BP-1 and BP-2. 
Hence, on an average, one would expect the leptons to be a little harder compared to the background leptons 
originating in the cascades of the top quarks. This can be seen from the left panel of figure [ 6 ] where the pt 
distributions of the harder lepton in the signal in both the benchmarks have extended tails compared to a similar 
lepton originating from the SM background. The signal distributions for BP-1 and BP-2 are similar because of 
similar values of top squark masses in the two benchmarks. The signal Xr distributions are different from the 
corresponding distribution for the SM background. This may be attributed to the much larger mass of the top 
squark (compared to rrit) and the presence of extra carriers of Xt, he., the bino-like MeV neutralino (y!)) and 
the active neutrino (i/ e ) emerging from top squark decays. It is evident from the right panel of figure [ 6 ] that 
the SM background can be effectively suppressed by applying a hard enough Xr cut, viz., Xt > 200 GeV. We 
have also constructed the dileptonic stransverse mass variable to see if the SM background can be suppressed 
further. The stransverse mass is a kinematic variable which is used to measure the ma sses of the pair produced 
semi-invisibly decaying heavy particles. The dileptonic stransverse mass is defined as jl4C| 


m T2 (PtiPt^t) = miny T= ^ + y2 


max| M t ( Pt ,pr)i M T(P t -Pt)] 


(38) 


where t\ and ti are two isolated leptons and pT is the total missing transverse momentum in the event and 
the transverse mass of the system Mt has its usual definition. Although, the stransverse mass is a standard 
variable used in the recent top squark search, however, the yield with Xt cut is more effective compared to the 
dileptonic M^ 2 cut for BP1 and BP2. 

• pp —> t 2 t 2 —> bbe + e~ (figure [3j) 

As has been pointed out earlier, t 2 could have a significant decay branching fraction to a bottom quark and an 
electron, which is a characteristic of such a scenario. Along with the enlarged phase space available to this decay 
mode, a moderately large coupling (~ A ' 133 = y^) does boost the decay rate. Naturally, we expect electrons 
(positrons) with high pt- In the absence of a genuine carrier of Xr in such a final state, low or at most a 
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PP ~> t\ti* —» bbt+t 
1= e,p 


13 TeV Background 

13 TeV BP1 - 

13 TeV BP2 - 


I 

> 

<D 


o 



1 ^ 


10" 1 


io- 2 


10" 3 



pf(GeV) | r r (GeV) 

_ 13 TeV Background 

0.2 pp-tttfi -*bM + t-g T 13 TeV BP1 - - 

i = e,p 13 TeV BP2 - 



50 100 150 200 250 

M"( GeV) 


FIG. 6: pr distributions of the harder lepton (left) and the^T distributions (right) for the background and the 
signals and dileptonic stransverse mass (in benchmark scenarios BP-1 and BP-2) in the 2 6-jet+2 lepton+^r final 

state arising from decays of t\ for the case p u < . 



FIG. 7: pr distributions of the harder electron (positron) from both background and signal (left) and $t 
distributions (right) (in benchmark scenarios BP-1 and BP-2) in the bbe + e~ final state arising from decays of a pair 

of tz-s for the case p, u < mj . 


moderate is expected from mis-measured momenta of the involved physics objects. The leptons are also 
expected to have uncorrelated momenta. Such events are rare in the SM. The left of panel figure [7] illustrates the 
hardest electron (positron) pt distribution in the scenario where t 2 decays to a bottom quark and an electron 
(positron). We impose a minimum pr cut of 200 GeV to reduce the SM background substantially. Since my 2 
is very similar for BP-1 and BP-2 and so is its kinematics for these two benchmark points, the distributions 
look very similar. In the right panel of figure [7] we present the MET distribution which arises in this case from 
mis-measurements of momenta of visible entities in the final state. As expected, the MET distributions peak at 
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PU GeV) 

_ 13 TeV Background 

0.2 - pp^t\t\ —► bb£ + £~Xr 13TeVBP3 - - 

0 _ 13 TeV BP4 - 



Zr(GeV) 


50 100 150 200 250 


M“( GeV) 


FIG. 8: pr distributions of the harder lepton (left) and the^T distributions (right) and dileptonic stransverse mass 
for the background and the signals (in benchmark scenarios BP-3 and BP-4) in the 2 6-jet +2 lepton-f^Er final state 

arising from decays of ti for the case p u > . 


small MET (« 25 GeV). 

Note that eventually, one should be able to reconstruct t 2 -s in the invariant mass spectra of appropriately chosen 
6 -jet-electron (positron) systems which would show peaks at rnj 2 . Clearly, the efficiency of reconstructing t 2 
would be limited by various detector effects and a close study of the kinematic distributions discussed above 
would surely be of crucial help. Nonetheless, it appears that the peaks cannot be missed and a reasonable 
estimation of mj would thus be possible. 


2. Case 2: p u > 2 


For p u > my , the top squarks decay mostly in a symmetric manner with t\ —> t \2 and 62 —> 6 e + as can be seen 
from table IIXI 

• pp — > t\t\ —> t \2 t\° —*■ 2 6 -jets + 2 leptons + Xt (figure [ 2 b]) 

For p u > mj , t\ decays mostly to a top quark and a bino-like neutralino (see table IVIIII) . The top quark 
would subsequently decay to a W boson and a 6 -jet via cascades. A pair of IP’s can then decay leptonically, 
semi-leptonically or hadronically. We confine ourselves to leptonic decays of IF-bosons for cleaner signals. The 
final state would then be comprised of 2 6-jets+2 leptons-bi^r- The pt distributions for the harder of the final 
state leptons are shown in the left panel of figure [5] The presence of an additional source of J?t and the heavier 
mass of t\ in the signal are behind harder distributions (see right panel of figure (JHI) when compared to 
the SM background. To optimise the signal significance, we have observed the dileptonic cut of 150 GeV 
works better compared to the fir cut of 200 GeV as used in Case 1. 

• pp —> 62^2 ~t bbe + e~ (figure [3]) 
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FIG. 9: Same as in figure [7] but for benchmark scenarios BP-3 and BP-4 and for the case p u > rrip 2 . 


In this case the overwhelmingly dominant decay mode is £2 —t 6 e + . As mentioned earlier, the emitted electron 
(positron) could have a very high pr as is evident from the left panel of figure [9] A strong pt cut (>200 GeV) on 
the electron can thus be easily afforded to suppress the SM background effectively. Similar to the case of figure 
0 the distributions of MET (of spurious origin) for the present case are presented in the right panel of figure 
[9] Again, the MET distributions peak at small values (~ 25 GeV), as expected and explained earlier. Again, 
possible reconstructions of t 2 -s in the invariant mass distributions of suitable pairs of 6 -jet-electron (positron) 
systems are on the cards. We would touch briefly on this issue later in this section. 

It is also important to note that ATLAS has performed a search for RPV stops in this channel [l4l| . Our 
analysis strategy is somewhat different from what they have chosen. ATLAS uses a cut on hadronic transverse 
momentum Hr > 1.1 TeV and requires the invariant masses of the 6 -lepton pairs to be within 20% of each 
other. We find that the simple cut on electron momentum suppresses the background equally well and should 
be robust even at high pileup conditions. 

Before going into the assessment of the signal significance, we mention below some issues of interest/importance 
pertaining to possible final states in these two cases. 

• For both p u < (section IVI AD and p u > rrij (section IVIBD . we have only looked into the 2 6-jet+2 
lepton+ / E l 7 ’ final state arising from t\ pair production. However, the first scenario is phenomenologically richer 
as it can yield multi-lepton signals with 4-6 leptons in the final states when t \-s and t 2 ~s decay via cascades 
involving the heavier neutralinos and charginos that in turn decay to SM Z bosons (see tables Hill and fVl) . Some 
corroborative analyses can take advantage of such inclusive final states comprising of 4 to 6 leptons along with 
6 -jets and MET. 

• Furthermore, such a possibility could help differentiate 62 from the two distinct scenarios considered in this work. 
For p u < my 2 , in addition to the 2 6 -jets + e + e - final state out of which a pair of £2 could be reconstructed, 
there would also be multi-lepton final states where leptons other than e + (e _ ) would appear. This is sharp 
contrast to the regime with p u > my 2 . 

• A final state like 2 6 -jets + e + e _ arising from the decays of t- 2 -s would be ideally free from any MET. However, as 
pointed out earlier, in reality, mis-measurements of various momenta may give rise to low to moderate amount 
of MET thus rendering the final state arising from a pair of t, 2 -s to be similar to that is obtained from t\ 
pair-production in a part of the phase space. This gives rise to some legitimate concern as to how efficiently 
the signature of t 2 -s could be deciphered, given the rates for such a final state originating in t\ pair production 
would be, in general, large thanks to smaller mass of t\. 

Such contaminations, however, can be avoided to a reasonable extent by imposing hard cuts on the minimum 
Pt of the leading electrons as guided by the lepton pt distributions in the left panels of figures [ 6 ] and 0 (for 
p u < rrit ) and figures [8] and 0 (for p u > m^). In addition, imposition of a cut on the maximum allowed MET 
could effectively restrict the contamination thus allowing for a more efficient reconstruction of 62 - By studying 
the MET distributions presented in the right panels of figures 0 and 0 we find an optimal value of this cut to 
be Et < 50 GeV that helps retain a healthy number of ‘signal’ events with low Et, a characteristic of such a 
final state originating in the decays of 62 - 
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FIG. 10: Invariant mass distributions for the appropriate pairs of 6 -jet-eletron (positron) systems (with low MET 
characteristic of A-parity violating decays of £2 to a bottom quark and an electron (positron)) for p u < m t ~ (left) 
and p u > rtif 2 (right). The distributions are obtained by imposing pt > 200 GeV for the leading electron and 

Xt < 50 GeV. 


In figure [10] we present the invariant mass distributions of appropriately chosen pairs of 6 -jet and an electron 
(positron). Guided by figures [7] and O a high pt threshold of 200 GeV for the leading electron is demanded 
along with requiring a Xr < 50 GeV to ensure that we mostly confine ourselves to the signal region. The left 
(right) panel of figure fTOl represent the case with p u < m t ~ (p u > mp 2 ). We find that in both cases clear peaks 
at m t ~ show up thus raising the hope that not only t 2 ~s could be discovered in this mode but also a reliable 
estimate of its mass would be possible. 

Before we close this subsection we like to mention that although we have only discussed two broad scenarios, i.e., 
Pu < m t ~ and p u > ro t ~, other intermediate situations are all a priori viable. However, the expectations under 
those scenarios could be substantiated in a straightforward manner from the two cases we present. For example, an 
increase in value of p u from that in Case 1 would result in suppression of the branching fractions to Higgsino-like 
neutralinos and charginos. With increasing /z u , at some point, these decay-modes (see table 3) would be closed 
for t\ and BR(ti —> t\ 2 ) = 1- At the same time, branching fractions to the Higgsino-like states for I 2 would also 
get suppressed before these decay-modes get completely closed as it happens in Case 2. A detail study of possible 
correlations among the event rates in various final states could, in principle, shed light on the relative value of p u 
with respect to mp and mp . However, this is beyond the scope of the present work. 


E. Signal significance and the reach 


The signal significance (a) is estimated using the expression [14: 

<7=^2[(S + £)ln(l + !)-S 


(39) 


which is appropriate for the situation with small number of events (in particular when the number of background 
events is less than 50). Equation (l3!?l) is based on likelihood-ratios and follows from the Poisson distribution. Here, 
S and B stand for the numbers of the signal and the background events, respectively after imposition of the set of 
optimal cuts discussed in section IVH Dl The A-factors for tt* are computed using Prospino2 (v2. 1 ) Gum. 

We now estimate the required integrated luminosities for a 5tr reach of t\ and £2 hr the four benchmark scenarios 
we consider. The final states we focus on are 2 6 -jets +2 leptons+Xr and bbe + e~, which stem from the decay of t\ 
and £ 2 , respectively. For the first case (see table m the dominant background comes from tt pair production which 
subsequently decays to the 2 6 -jet +2 lepton+Xr final state. An appropriate A'-factor of ~ 1.6 is used to derive the 
NLO cross sections from the LO ones for f-pair production. 

We note that a 5<r signal significance can be achieved for BP-1 and BP-2, with an integrated luminosity around 100 
fb” 1 . To achieve a similar significance for BP-3 and BP-4, one has to wait for a much higher accumulated luminosity, 
for example, 500 fb -1 and 3000 fb -1 , respectively at yfs = 13 TeV. ~ 

In table m we present the required luminosities for a 5<r reach of £2 for the four benchmark points. The dominant 
SM background comes from bbZ/ 7 * production followed by Z/'y* giving rise to e + e _ pairs. This can be efficiently 
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pp — > titj —¥ 2 6-jets + 2 leptons +J?t 

BP1 

BP2 

BP3 

BP4 

a(pp^titl) (fb) 

428.9 

463.4 

193.0 

73.6 

Cut acceptance for signal 
(J?r > 200 GeV for BP1 and BP2) 
(M|4 > 150 GeV for BP3 and BP4) 

1.5 x 10” 2 

1.6xl0” 2 

4.4x 10” 3 

4.4 x 10” 3 

Required £ (fb^ 1 ) for 5cr significance 

256.0 

316 

2350.0 

3000 (3a) 


TABLE X: Required values of integrated luminosities (£) to obtain a 5cr significance in the final state at yfs = 13 
TeV. The most important SM background arising from tt pair production is normalised to a cross section of ~ 816 
pb obtained at the NLO+NNLL level (see section lYII Cl) . The cut acceptance for the background is 2.3xl(V 4 . A 

flat 6-tagging efficiency of 60% is used. 


pp t 2 t 2 2 6-jets + e + e 

BP1 

BP2 

BP3 

BP4 

a(pp -*• t 2 tl) (fb) 

7.83 

8.69 

8.45 

8.48 

Cut acceptance for signal 
(p T > 200 GeV, J? T < 50 GeV) 

1.96xl0~ 2 

5.26xl0 -2 

1.9xl0 _1 

1.9xl0 _1 

Required £ (fb^ 1 ) for 5cr significance 

501.23 

63.92 

6.85 

6.61 


TABLE XI: Required values of integrated luminosities (£) to obtain a 5<r significance in the 2 6-jets + e + e~ final 
state at y/s = 13 TeV. The SM background (see section IVlI Cl) at NLO is found to be 9.43 pb (see section IVlI Cl) . 
The cut acceptance for the background is 1.25xl0 -4 . A flat 6-tagging efficiency of 60% is used. 



FIG. 11: Density plot reflecting the reach for top squark masses via 2 6-jet+2 lepton+Ffr final state in BP-2 in the 
my | -£ plane (left) and via 2 6-jet+ e + e _ final state in BP-4 in the mj r -L plane (right). The thick black curves are 
contours of 5cr significance above and on the left of which the masses can be explored with > 5 a significance. The 
cuts are kept fixed at values mentioned in the text for the respective cases. 


suppressed by using an on-shell Z-veto for the e + e _ pairs, as discussed in section I VII Cl Thus, as can be seen from 
this table, a 5cr significance can be obtained with an integrated luminosity as low as < 10 fb” 1 for the benchmark 
scenarios BP-3 and BP-4 with y/s = 13 TeV. In addition, we also study the Hp distribution, i.e., the scalar sum of the 
Pt of the e + e” pair and the reconstructed 6-jets and the improvements are marginal. The wildly varying integrated 
luminosities across the benchmark points are the artifact of varying branching fractions that are instrumental, as has 
been pointed out in section ED 

Figure QT] summarises the mass-reach for the two top squarks with varying accumulated integrated luminosities (or, 
in other words, luminosity required to probe a certain top squark mass) at the 13 TeV LHC. The left panel illustrates 
the case for t\ in the final state 26+2 lepton in BP-2 while the right one does the same for 1 2 via 26+e + e” final 
state in BP-4. As indicated by tables IXl and IXI1 figure HT1 also reveals that t 2 has a significantly better reach compared 
to t\ with the final states under consideration. This may lead to a tantalising possibility of discovering t 2 of such a 
scenario much earlier than t\ and the former could guide us to find the latter. We observe that at the 13 TeV LHC 
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FIG. 12: pt distributions of the harder lepton (left) and the t distributions (right) for the background and the 
signal in the 26-jet+2 lepton + final state arising from the decays of a pair of t i-s in the ‘stealth’ top squark 

scenario. 


and with the mass-reaches for t\ and t-2 are around 575 (750) GeV and 1.2 (1.4) TeV respectively, with an integrated 
luminosity of 300 (3000) fb -1 . 

VIII. THE ‘STEALTH’ TOP SQUARK SCENARIO 

The SUSY model under consideration, with super-light carriers of MET like Xi = and an MeV neutralino LSP 
(X 2 ) can easily conceive a rather low mass top squark lying right in the so-called ‘stealth’ window of 197 GeV < < 

205 GeV [ 73 I. As discussed in section IV Cl the experimental lower bound on mp L is more stringent considering its 
decay modes. Hence we choose t r to be the lightest top squark (U). A benchmark point can be obtained by choosing 
( to u) 33 = —2.5 x 10 4 GeV 2 . This results in mj ~ 200 GeV. Such a light top squark cannot provide enough correction 
to the Higgs mass. Hence we choose a relatively large value of As (=1.28) so that the radiatively generated additional 
quartic contributions could lift the Higgs boson mass to the observed range. All other parameters are fixed at the 
values mentioned in BP-1 (see table HT1) . Note that the additional tree level contribution proportional to the neutrino 
Yukawa coupling ‘/’ remains small (even for its order one value) because of large values of tan/3 that we require. As 
a result, t\ mostly decays to t \2 and ti —> tv e with ~ 85% and ~ 15% branching fractions, respectively. The possible 
final state topologies are exactly the same as those result from top quark pair production. 

We again analyse the final state with 2 6-jets+ 2 leptons-hiSr- We checked that the distributions of various kinematic 
observables look very similar for the signal and the ti background, which is something literally expected of a ‘stealth’ 
top squark and what makes it so elusive. In figure [12] we present the pt distribution of the harder lepton (left panel) 
and the distribution (right panel) which clearly demonstrate how similar the behaviors of the SM background 
and the signal could get. In this context, techniques to e xploi t differences in spin-correlations inherent to ti and 
tt* syste ms /82| . use of various transverse mass variabl es Il43l including the one like ttit 2 in the dileptonic decay 
channel |l44j . incorporating a new variable like ‘topness ’ m using asymmetric decays of the top squarks have been 
proposed to study the ‘stealth’ top squark regime in search for an improved sensitivity. Clearly the issue demands 
dedicated addressal which is beyond the scope of the present discussion. 


IX. SUMMARY AND CONCLUSIONS 

We study a £7(1) r —lepton number model augmented with a right handed neutrino superfield. The i?-charges are 
identified with the lepton numbers in such a way that the sneutrinos acquire large vev-s. Such large vev-s for sneutrinos 
are not prohibited since, in such a scenario, the same are not constrained by the Majorana masses of the neutrinos. In 
this paper, we choose to work in a basis in which only the electron type sneutrino acquires a nonzero vev whereas the 
vev-s of the other two sneutrinos are rotated away. This simple extension with a right handed neutrino superfield is 
rather interesting in the sense that the Higgs boson mass gets a tree level contribution which can be substantial in the 
low tan/3 regime and for order one neutrino Yukawa coupling ‘/’. Also present are the large one loop contribution to 
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the Higgs boson mass arising from new couplings in the theory. Thus, one can easily accommodate a scenario where 
both top squarks are light. In addition, a very light bino-likc neutralino comes out naturally in this scenario along 
with an active neutrino endowed with an appropriately small Majorana mass. Therefore, rich and interesting collider 
signatures are expected in such a scenario. The signature of a top squark decaying to a top quark and a neutralino 
and/or a neutrino is similar to the top quark pair production in the SM. Under favorable circumstances, top squark 
could also decay to a bottom quark and a chargino leading to a similar final state containing 2 6 -jets and 2 leptons 
along with MET. In this mode, top squark mass of around 575 (750) GeV can be probed with 300 (3000) fb -1 of 
integrated luminosity. 

Furthermore, in an 17-parity violating scenario such as ours, the cliarginos mix with the electron. The decay width 
of the top squark to a bottom quark and an electron (positron) is enhanced because of the enhanced coupling A'^g as 
well as an unsuppressed phase space. Hence we study in detail the final state with 2 6 -jets accompanied by an e + e - 
pair arising from such a dominant decay. We show that even when the top squark is heavy (jnj 2 ^ 1.2(1.4) TeV), this 
particular channel could deliver a large signal significance with 300 (3000) fb _1 of integrated luminosity. 

In the model discussed in this work, only i 2 (~ ti,) decays to 6e + . This is an artifact of no L-R mixing in the top 
squark sector. Also, note that the scenario prohibits G decaying to 6/i + or 6r + . The final state arising from both 
top squarks decaying to be + {be~) mode could carry MET which can only be of spurious origin (mis-measurements 
of various visible momenta, defects in the detector, etc.) and hence is characteristically small. This feature can be 
used to establish such a model and differentiate it from other competing ones. Due to a relatively clean final state 
and hence, a possibility to reconstruct the heavier top squark mass reasonably efficiently, such a state could be within 
an easier reach of the current LHC run when compared to its lighter peer. Such a scenario thus, gives rise to an 
interesting possibility that t >2 can be found much earlier than t\ at the LHC and could carry a reliable hint as to where 
exactly to look for the latter. The signal region for t\ is attributed with a much larger MET as is usual in searches 
for new heavy states in scenarios with a stable charge- and color-neutral particle(s). This is in sharp contrast with 
the case of G in such a scenario. 

Although the analyses in this work are presented in terms of two broad scenarios, viz., /i u < m t ~ and fi u > m f ~, 
it is pointed out that the signatures discussed are robust under intermediate situations except for some obvious 
quantitative issues getting in. Simultaneous searches in various channels described in this work are expected to shed 
light on the detailed aspect of the spectrum and the involved new couplings of such a scenario. 

Finally, we have demonstrated how the ‘stealth’ top squark can appear in our model naturally. However, probing 
such a window, 197 GeV < mj- < 205 GeV [78] needs dedicated analysis which is an active area of research on its 
own merit. Overall, characteristic signatures for these light top squark states at the LHC even have the potential to 
discriminate between competing scenarios that may give rise to such a light pair of top squarks. In addition, such 
issues and projections are not readily available for 13/14 TeV run. Therefore, it is important to study all these issues 
at the dawn of 13 TeV run of the LHC. 
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